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An Analysis Of Optical Pick-up In SCUBA Data 
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ABSTRACT 

The Sub-millimetre Common User Bolometer Array (SCUBA) at the JCMT employs a chop- 
ping and nodding observation technique to remove variations in the atmospheric signal and 
improve the long term stability of the instrument. In order to understand systematic effects in 
SCUBA data, we have analysed single-nod time streams from across SCUBA's lifetime, and 
present an analysis of a ubiquitous optical pick-up signal connected to the pointing of the sec- 
ondary mirror. This pick-up is usually removed to a high level by subtracting data from two 
nod positions and is therefore not obviously present in most reduced SCUBA data. However, 
if the nod cancellation is not perfect this pick-up can swamp astronomical signals. We discuss 
various methods which have been suggested to account for this type of imperfect cancellation, 
and also examine the impact of this pick-up on observations with future bolometric cameras 
at the JCMT. 
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1 INTRODUCTION 

In recent years, millimetre (mm) and sub-millimetre (sub-mm) 
bolometric cameras have become the focus of much effort. These 
instruments directly detect thermal emission associated with star 
formation and have proven useful in studying previously hid- 



shifts (e.g. 


Lilly et alll999llDunne & Ealesll200 ill Johnstone et alJ 


2001, Ivisonetal. 


2002, 


Scott et alJ 


2002, 


Chapman et al. 2003). 



The S ubmillimetre Common User Bolometer Array (SCUBA; 
iHolland et alJ 1 19991) at the James Clerk Maxwell Telescope 
(JCMT), with its large aperture and excellent site, has been par- 
ticularly productive in this field. 

Even exploiting a combination of high, dry sites and 'win- 
dows' in atmospheric opacity, emission due to the atmosphere 
and telescope's surroundings is very large compared to astrophys- 
ical signals ( Archibald et al 1 120021) . Moreover, it is highly vari- 
able on time scales slower than a few Hz. To counteract this, 
sub-mm telescopes can employ a variety of observing strategies 
to scan the sky rapidly and take differences to remove the bulk 
of the atmospheric emission. These techniques put astrophysi- 
cal signals into an audio frequency band above the atmospheric 
1/f knee. Examples of specific observing s trategies are those em- 
ployed by SHARC II iDowell et aill2003l). which sca ns the sky 
in a lissajous pattern, and BLAST l Tucker et al. 2004) and ACT 
iKosowskv) 120041) . which will scan the entire telescope in az- 
imuth only. SCUBA uses a rapidly moving secondary mirror to 
chop the beam rapidly and nods the whole telescope more slowly 
( Holla nd etalll999l) while SCUBA2 (which will also be deployed 



at the JCM T.lHolland et alj|2003h intends to emp loy the DREAM 
algorithm iLe Poole & van Someren G reve 1998). In all of these 
approaches the observation and data reduction processes which 
cancel atmospheric signals also filter astrophysical signals at large 
angular scales. 

Previous efforts to measure the Sunyaev-Zel'dovich (SZ) ef- 
fect on the JCMT have failed due to early removal of large angular 
scale signals during data reduction. Therefo re, as part of a n effort to 
measure the SZ effect with SCUBA IZemcov et aiT2 003). we have 
taken the unorthodox step of considering SCUBA data before sub- 
tracting data from adjacent nods. When we do this we find a repet- 
itive signal which is orders of magnitude larger than typical astro- 
nomical and even atmospheric signals. Fig.Qshows an example of 
this signal. We will show that this signal is associated with changes 
in the orientation of the secondary mirror and presumably arises 
from altering the amount of thermal emission from the telescope 
which is received by the bolometer array. Although this signal is 
present in all SCUBA data we have examined, it generally seems to 
be rem oved effectively in the standard SURF Jjenness & Lig htfoot 
1998) reduction pipeline. However, a period during which nod sub- 
traction does not effectively cancel this pick-up has already been 
identified by at least one author <Borvs et all2004l) . 

This signal is synchronized with motion of the secondary mir- 
ror, and therefore occurs in the data at a set of audio frequencies 
overlapping those of astronomical signals. Its presence imposes 
severe constraints on viable observation strategies and on the re- 
quired precision of optical control at both the JCMT and similar 
telescopes. In this paper we present measurements and analysis of 
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Figure 1. A typical reconstructed SCUBA data series for a 120 arcsecond 
chop throw. This particular data set is for bolometer G13 from observation 
28 on 1998 June 30. The upper trace shows the time ordered data collected 
during four consecutive 64 point jiggle patterns. Each sample is 1.024 s of 
data averaged over the chop but there are occasional small gaps in the data 
series associated with time spent nodding the telescope not shown here. Ev- 
ery bolometer in the array shows a similar pattern which is large, stable, 
and occurs at the jiggle frequency. In the usual data reduction pipeline suc- 
cessive nod patterns are subtracted and this cancels the repetitive pattern 
almost completely. This is shown as the blue interrupted trace which repre- 
sents Al - Bl, A2 - B2, A3 - B3, and A4 - B4 for the four cycles. 
There are therefore half as many data points in this double-differenced set 
as in the raw data, so this trace is plotted only in the A nod position. The 
lower red (continuous) curve is the residual after subtracting the polynomial 
model of mirror position discussed below. It is nearly as clean as the blue 
data set. 

this signal with the aim of learning what these constraints are and 
determining the origin of this signal. 
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Figure 2. A schematic of the SCUBA jiggle pattern. It is controlled by 
motion of the secondary mirror, and allows the SCUBA array to fully sam- 
ple the sky. The central dashed circle shows the full width half maximum 
(FWHM) beam size associated with the (circled) central triangular point. 
The peripheral dashed circles show the FWHM of adjacent bolometers in 
the SCUBA array; SCUBA under samples the sky at any instant. The sec- 
ondary mirror chops the beam on the sky at ~ 8 Hz between positions 
separated by a user-defined angle, dwelling for ~ 1 s with one end of the 
chop on each of the positions shown above. After completing the 16 posi- 
tions in one quadrant of the pattern the telescope primary mirror nods and 
the other end of the chop is put onto each of the 16 positions in the same 
quadrant. Then the second, third, and fourth quadrants of the jiggle pattern 
are observed in the same way. If the primary mirror nod positions are called 
A and B and the jiggle pattern quadrants are numbered 1 through 4 as in the 
figure, a 64 point measurement pattern consists of the following sequence 
of 16 s measurements: Al Bl B2 A2 A3 B3 B4 AA. This process requires 
approximately 128 s to complete, plus some overhead associated with mov- 
ing the primary mirror. Notice that the nod frequency of the primary mir- 
ror is ~ 1/32 Hz while secondary mirror positions occur at harmonics of 
~ 1/128 Hz 



2 OBSERVATION OF OPTICAL PICK-UP 
2.1 SCUBA observation strategies 

iHolland et all \ 19991 and lArchibald et alJ l2002h provide excellent 
and extensive reviews of standard SCUBA operations and system- 
atic effects. We review here just enough of the SCUBA observation 
modes to allow understanding of the pattern shown in Fig.Q 

The JCMT's beam is chopped on the sky by nutating the sec- 
ondary mirror at about 7.8 Hz in a nearly square wave pattern. The 
angular separation of the end points of this chop, called the chop 
throw, is chosen by the observer. The chop is symmetric about the 
optic axis of the telescope, which is to say that the two endpoints of 
the chop pattern are each displaced from the optic axis by half the 
chop throw. The measured difference in intensity between the two 
endpoints for every bolometer in both arrays is reported at approx- 
imately 1 Hz. This signal is largely free of the uniform component 
of atmospheric emission. We will call this differential data the raw 
time series. 

The entire telescope is reoriented periodically to put one end 



of the chop pattern and then the other on to the nominal source 
position. This is called nodding the telescope. The difference in 
intensity between the two nod positions forms a triple difference 
on the sky: the brightness at the nominal source position minus that 
at two symmetrically located off-positions. This signal is largely 
free of the effects of gradients in atmospheric emission. This data 
set is called the de-nodded or double-differenced data, and is also 
shown in Fig.Q 

The SCUBA array under fills the focal plane, so the secondary 
mirror undergoes a slow, small amplitude walk in addition to the 
movements described above so that over time the sky is fully sam- 
pled. Fig. [2] shows the standard 64 point pattern, called a jiggle 
pattern. The 64 points are split into four quadrants of 16 pointings 
each and the telescope nods after observing each quadrant. If we 
label these quadrants 1 through 4, and label the two nod positions 
as A, and B, the order of data collection is 

Al Bl B2 A2 A3 S3 Bl Al Al . . . 

The pattern repeats every (16 Positions) x (2 Nods) x 
(4 Quadrants) = 128 samples, requiring slightly longer than 128 
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seconds to complete one cycle. The sample series also has gaps in 
time because the primary mirror requires time to move. The order 
of observing nod A and nod B alternates by quadrant to minimize 
primary mirror motion. The peak to peak jiggle distance is about 30 
arcsec, and the average jiggle distance from the middle of the pat- 
tern is about 9 arcsec. Very similar, albeit simpler, motions occur in 
SCUBA's 16 point jiggle pattern and 9 point photometry modes. 

Although the SCUBA bolometers are sampled at 128 Hz, data 
are reported as one second averages over the chop; undifferenced 
data are not saved. In standard operation, two time streams are re- 
ported, one for each nod position. Normally, the nod time streams 
are subtracted from each other (Al — Bl etc.), thereby removing 
a large component of the atmospheric emission and instrumental 
variation. This double-difference time stream, which consists of 
the target beam data minus half of each of the two chop positions' 
data, is the standard SCUBA data product. 

2.2 The jiggle frequency signal 

The raw time series in Fig.Qshows a clear periodicity at 128 sam- 
ples, indicating emission is modulated by the 128 point jiggle-and- 
nod pattern. The peak to peak amplitude of the pick-up signal is 
about 10 mV, which is about 10 4 times as large as those arising 
from the weakest astronomical sources SCUBA is capable of re- 
solving. In order to provide a consistent, if approximate, guide by 
which to compare these signals to astrophysical sources and ther- 
mal effects we have adopted a gain of 210 Jy/V throughout this pa- 
per. This is consistent with the mean value observed in pre-upgrade 
SCUBA data Jjenness et all20 02). With this gain, the amplitude of 
the pick-up is about 2 Jy per beam peak to peak. 

In the de-nodded data set there are 16 points in the set 
Al — Bl, 16 in the set A2 — B2, and so on. These points are 
also plotted in Fig.Q The amplitude of the pick-up signal has been 
dramatically reduced and these signals are consistent with a com- 
bination of detector and atmospheric noise. We conclude that the 
source of the pick-up is not modulated by nodding the telescope, 
and therefore that the pick-up signal is associated with motions of 
the secondary mirror only. 

To check this conclusion, and to quantify the relation, we fit 
the raw time series to a polynomial in a secondary mirror orienta- 
tion using elevation, h, and azimuth, A, coordinates as tracked by 
bolometer H7's jiggle offsets, which are nominally zero on the op- 
tical axis of the telescope. This basis forms the natural coordinate 
system for exploring a signal which arises within or around the 
telescope. We perform a fit independently in h and A, using each 
angle and its derivatives. It is found that the fit does not require 
cross terms of A and h to describe the data. Our fitting function is 
therefore 

9 9 9 

n—l n—l n—l 

9 9 9 

J2S n h n + ^2e n h n + ^2Cnh n (1) 

n—l n—l n—l 

where (A(t), h(t)) is the secondary mirror pointing at time t mea- 
sured in arcsec, dots denote time derivatives, and the constants a n , 
/3n,fn,S n , e„ and fn are determined from the fit. In practice, these 
fits are performed independently to A and then h for each bolome- 
ter. The number of terms in the fit is arbitrary; we have chosen 
n ^ 9 as the ninth order terms are consistent with zero. If the 
signal is optical in origin, we expect the n — l term, which corre- 
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Figure 3. Power spectrum magnitudes (PSM, the modulus of the Fourier 
transform) of the data shown in Fig.Q The Fourier transforms have been 
taken with respect to sample number as the precise mapping between sam- 
ple number and time is complicated due to small gaps in data collection 
during normal SCUBA operation. As a guide, 1 sample ~ 1 Hz. The black 
curve shows the PSM of the raw data set. Note that there is no power at 1/32 
sample -1 , confirming that the pick-up signal is associated with secondary 
mirror position and not with motion of the primary mirror. The red curve is 
the PSM of the residual signal after removal of the polynomial model dis- 
cussed below. This curve contains no power in any of the spikes in the raw 
data spectrum but faithfully reproduces the PSM everywhere else. The blue 
curve is the PSM of the double differenced data. Above 1/64 sample -1 the 
PSMs all have comparable power, but at low frequency double differencing 
produces 'whitened' data with about half of the variance of the model sub- 
traction method. 



sponds to position, to dominate, while if it is due to pick-up from 
motors or electrical equipment, higher order terms, corresponding 
to velocities and accelerations, should dominate. Any pick-up as- 
sociated with the nod does not affect the time series at the jiggle 
frequencies, but rather adds a variable DC offset between the two 
nod positions. When Eq. Qis fit, the linear terms ai and 8i are 
at least an order of magnitude larger than any other term. This fit 
is therefore consistent with optical pick-up alone; the higher or- 
der terms (n > 1) are included in order to demonstrate that the 
pick-up pattern can be completely described in terms of the po- 
sition of the secondary mirror. Further, we find that the pick-up 
pattern is not constant but varies smoothly with bolometer position 
(as shown below), which strongly suggests the pick-up is optical 
and not electrical in origin. Residuals from this fit, V(t) — V(t), 
are shown in Fig. Q Just like the de-nodded data, these residuals 
also show nearly complete elimination of the pick-up signal. In a 
single data file this cancellation is effective to about one part in 
100; unfortunately, to observe weak astrophysical signals in multi- 
ple integrations, we require cancellation to better than one part in 
3000. 

Fourier transforms of the raw data, the de-nodded data, and 
the residuals from Eq. Qare shown in Fig. [3] Caution is required 
when interpreting the transform of the de-nodded data since that 
data set is sampled differently in time than the other two sets. The 
raw data shows spikes at harmonics of 1/(128 samples). Notice 
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that the fourth harmonic has no amplitude above the continuum 
level. This is because the transitions from Al to Bl, A2 to B2, 
etc. occur at the fourth harmonic and nodding does not alter the 
secondary mirror pattern of motion. All of the harmonic spikes are 
missing from the power spectra of either the residuals or the de- 
nodded data. Except at the spikes, the fit residuals match the raw 
data with high fidelity. Notice that at very low frequencies the de- 
nodded data is a bit quieter than the residuals are. Presumably noise 
and large angular scale astronomical signals which occur at low 
frequencies are removed from the data by nod subtraction. 



3 CHARACTERIZATION OF THE PICK-UP SIGNAL 

In order to understand the severity and prevalence of this signal, 
data were obtained from both the SCUBA archive 1 , and a dedicated 
set of observations with very small chop throw in 2003 June. These 
data were chosen to reflect a sampling of the types of mapping ob- 
servation commonly performed with SCUBA. We have examined 
data taken in both 16 point and 64 point jiggle patterns, as well 
as SCUBA's photometry mode which uses a small 9 point grid. A 
variety of chop throw angles between 3 and 180 arcsec have been 
chosen. In order to reduce possible contamination by astronomical 
sources, the data sets were selected based on low estimated target 
flux. An attempt was made to use only fields containing sources 
with fluxes less than lOmJy, although this criterion was not met 
in the very small chop throw (less than 15 arcsec) data. A wide 
range of chop throws are chosen to give a wide sampling of pos- 
sible SCUBA jiggle map configurations. Also, data straddling the 
1999 October SCUBA refit and upgrade are considered. All of the 
data were taken on photometric nights, 0.04 < rcso < 0.09. In 
all, over 200 data files of varying lengths spanning nearly all of 
SCUBA's operational lifetime have been considered in this analy- 
sis. 

The data sets considered here have been preprocessed using 
SURF 2 . Programs are called to flat field, de-spike and extinction 
correct both the double-differenced and single nod time streams; 
we direct the reader to SURF's documentation for the details of how 
these procedures are implemented. All of the data discussed here 
are left at this stage of processing; advanced atmospheric removal 
has not been implemented. 

3.1 Shape of the pattern 

For each bolometer in a given observation, the vector 

g = a lv A + 8 lv h, (2) 

where the subscript r\ denotes the different bolometers, is closely 
related to the gradient of V in azimuth-elevation coordinates. Be- 
cause the linear terms dominate, g captures almost all of the power 
in V. The units of g are V arcsec -1 . When the vector field g is 
plotted against the location of each bolometer in the focal plane, a 
clear pattern emerges (Figs.|4]and|5j. There is typically a null in 
the patterns and the magnitude of g grows with separation from the 
location of the null in the focal plane. 

In order to test whether this signal is only present in the long 

1 Courtesy the Canadian Astronomy Data Centre, operated by the Do- 
minion Astrophysical Observatory for the National Research Council of 
Canada's Herzberg Institute of Astrophysics. 

2 For a detailed description of SURF, see http : / /www . starlink . 
rl.ac.uk/star/ docs/ sun216 .htx/sun216 . html . 



wavelength SCUBA array, the analysis discussed above is also per- 
formed on the short wavelength array (Fig.|4j. It is found that the 
patterns are remarkably similar between arrays, and in particular 
the null in the pattern occurs at the same position in the two focal 
planes during any given observation. This implies that this signal 
originates somewhere in SCUBA's optical path. 

The pattern in Fig.|4]is a derivative of the optical loading on 
each bolometer. There appears to be both a divergence and a curl 
in the pattern, perhaps because Eq. Q does not employ true par- 
tial derivatives, or because the secondary mirror's orientation and a 
bolometer's location in the focal plane are not identically related. 
In any case, the location of the null in the pattern corresponds to 
the location of an extremum in the optical loading of the bolometer 
array. It is striking that in our experience this extremum does not 
occur at the array centre, which we presume is the symmetry axis 
of the telescope. However, since the telescope is very difficult to 
align to arcsec accuracy 3 , it is not unreasonable to expect the sym- 
metry axis to deviate from the array centre by a few tens of arcsec 
(P. Friberg, private communication). 

As a routine part of our Sunyaev-Zel'dovich effect measure- 
ments, an aluminum reflector was placed at the cryostat window 
and observations with large chop throw were performed. This elim- 
inates all optical inputs to the bolometer arrays but leaves electrical 
pick-up from the various motors and acoustic pick-up intact. The 
resulting pattern is shown in Fig. [5] The average amplitude of g is 
more than two orders of magnitude smaller than without the reflec- 
tor in place. The pick-up signal must be optical in origin. 

The correlated g patterns occur in every SCUBA data set, 
and gross features such as the magnitude of the vectors are con- 
stant. However, details of each 'vector field', such as the centroid, 
curl, divergence, etc., vary slowly in time. For instance, positional 
changes in the vector field centroid of tens of arcsec can be tracked 
over an evening of observation. Furthermore, duplicate observa- 
tions of the same field performed months apart result in patterns 
which are quite dissimilar. We believe that this is because a given 
field's altitude and azimuth have changed over the intervening time, 
even though in sky-based coordinates the field's position is con- 
stant. We have had difficulty locating data taken with precisely the 
same altitude, azimuth and chop throw separated by a suitable inter- 
val of time, but surmise that such data would produce very similar 
patterns if this is an optical pick-up related to ambient radiation in 
the telescope dome and its surroundings. 

3.2 The { g | ) - chop throw distance relation 

Careful analysis of Fig.|4]shows that the mean magnitude of the g 
vectors, (|g|), is correlated with chop size. The quantity (|g|) and 
its standard deviation as a function of chop throw are plotted in 
Fig.|6|for a sampling of data files with different chop throws. 

It is observed that the chop size and (|g|) are correlated. While 
the characteristics of this correlation at modest or large chops are 
relatively well defined, the signal's behavior at very small chops is 
not so clear. In an attempt to clarify whether (|g|) tends to zero at 
small chop throws, power laws with an offset, 7 

(|g|) = a/+ 7 . (3) 

are fit to data collected with chop throws <j> less than 30 arcsec. The 
parameters a and /3 are determined from the fit. The results from 

3 This is performed either by using an alignment laser or, post-2002, by 
examining the noise properties of the atmosphere. 
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Figure 4. The set of g for a variety of chop throws (quoted in the upper corner of each plot). The top two plots are comparisons of the 450/850 (im arrays. 
The plot on the upper left shows the long wavelength array g from observation 67 on 1999 November 10, and the plot on the upper right shows its short 
wavelength counterpart. The g are shown beginning at each bolometer location ('+' symbols). Axis units apply only to the bolometer locations, representative 
scales can be found in the lower left of each plot. Bolometer 1 is marked with a triangle symbol and a scaling for g is shown in the bottom left of each plot 
for reference. The patterns are remarkably similar for both arrays, suggesting that the jiggle frequency signal originates in SCUBA's optical path rather after 
the dichroic beam-splitter i Holla nd et alll999h or in the electronics. The bottom two plots are for a 68 arcsec chop (observation 119 from 2003 June 14; left) 
and an 8 arcsec chop (observation 143 from 2003 June 14; right). In 2003, bolometer G9 (number 7) was not operational; its data from that period are not 
included. These lower plots should be compared with the 180 arcsec chop throw plot; on a gross scale, the patterns seem to be similar over time. However, 
smaller chops result in smaller pick-up signals, and the exact pattern changes with each observation. 



two of these fits, one with an offset of 7 = V arcsec - and the 
other with a voltage offset of 7 = 2.5 x 10 -5 V arcsec -1 (this value 
being arbitrarily chosen for illustrative purposes) are listed in Table 
1 and plotted in Fig.|6| The x 2 statistics of these fits are essentially 
the same. Although the two functions differ in principle, for chop 
throws of one arcsec the range allowed for the pick-up amplitude is 
small, from 5 to 7 mjy per arcsec. Even if the pick-up pattern tends 
to zero at (j> — 0, coherent mis-pointing during a jiggle observation 
by a few per cent of the JCMT's FWHM results in a spurious signal 
comparable to extra-galactic signals of interest. 



Table 1. Power law fit parameters. 



Fit 


7 


a 




P 




(V arcsec -1 ) 


(V arcsec - 


x ) 


(arcsec -1 ) 


'No Offset' 


0.0 


2.07 X 10 


-5 


0.586 


'0.5 mjy arcsec -1 ' 


2.5 x 10 -5 


7.35 x 10 


-6 


0.845 



Interestingly, the power laws hold even for the large chop data, 
although using data at large chop is a poor way to constrain the 
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Figure 5. SCUBA pattern with a reflector blocking the dewar aperture (ob- 
servation 26 from 2002 October 10), with symbols as in Figure 3. The g are 
consistent with being random, supporting the claim that the jiggle pick-up 
originates outside of the SCUBA cryostat. 




Reflector Limit 
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Figure 6. Jiggle frequency signal magnitude as a function of chop throw. 
The experimental points and their standard errors are shown as diamonds 
with bars. The lines show power law fits to the smallest four chop throw data 
points. The best-fitting power law with 7 = V arcsec -1 offset is shown 
as the solid line. Also, a power power law with 7 = 2.5x10 V arcsec 
(~ 0.5 mjy arcsec ) offset is shown as the dot-dash line. The dashed line 
and arrow at bottom indicate a lower limit to this signal as determined from 
the reflector data. Because the reflector data's g are uncorrected, (|g|) is 
an overestimate of the effect in this case, and we plot the limit as half of that 
calculated. Clearly, the jiggle frequency signal has a substantial amplitude 
even at small chops. It is very difficult to ascertain the behavior of this signal 
as the chop throw approaches arcsec. 



characteristics of the signal at small chop (this is the motivation for 
only using the smallest chop data in the fit). While the presence of 
an offset cannot be ruled out by these data alone, we suspect that 
the relation is essentially linear and tends to zero at zero chop. 



4 DISCUSSION 

Finding an optical pick-up signal which depends upon the orienta- 
tion of the secondary mirror is not a surprise. Moving the secondary 
mirror alters the electromagnetic cavity formed by the primary and 
secondary mirrors of the telescope, and also alters the fraction of 
the beam which intercepts the telescope structure or misses the pri- 
mary mirror. As has been pointed out, the pattern in Fig. [4] is a 
derivative of total optical loading on the bolometers. It is quite strik- 
ing that the observed pattern sometimes corresponds to a minimum 
loading near to the optical axis, and sometimes to a maximum (as 
in panels 1 and 3 of Fig.|4j- 

To compare the observed signals to expectations based on a 
physical understanding of the telescope first requires conversion 
of the observations from Jy per beam to an equivalent Rayleigh- 
Jeans (RJ) temperature. The RJ brightness of a source at 850 /im is 
3.8 x l(n 17 W m~ 2 Hz" 1 str -1 K _ \ and the solid angle of the 



JCMT's primary beam is fl 
so 



1.13 x FWHM = 5.8 x 10" a str, 



= 5.8 x 1CT 9 str 

x 3.8 x l(T ir W m~ 2 Hz -1 K _1 str -1 = 

22.0 JyKT 1 . (4) 



Ideally, this conversion fac tor would be ca lculated as an integral 
over the SCUBA passband I Borvs. Chapma n & Scottll999T) . which 
alters the numerical value by a few percent. 

The 2.5 Jy peak-to-peak signal in Fig. Q corresponds to a 
115 mK RJ signal. A more useful number for understanding the 
emission mechanism is that the chop throw dependent emissivity 
coefficient (|g|) for a typical chop throw of 60 arc seconds is 



,|geo" 



= 50 mjy arcsec 1 = 2.2 mK arcsec 



(5) 



Given that the telescope temperature is typically 270 K, this cor- 
responds to a variation in emissivity of 10 -5 per arcsec, averaged 
over the array. 

Several mechanisms could plausibly give emissivity variations 
which are this high, and we do not know the SCUBA optical illu- 
mination pattern well enough to make a definitive distinction. It is 
certainly possible that several processes are comparable in magni- 
tude. The aluminum reflective surfaces are a few per cent emissive. 
Their apparent emissivity drops as the mirrors are tilted with re- 
spect to each other giving rise to a maximum emission on the opti- 
cal axis and a variation which is to leading order quadratic in angle. 
This mechanism will not produce a minimum on axis. However, the 
measured signal is the difference between the gradients at two chop 
positions, and thus may have either sign depending on the order of 
the difference. 

The SCUBA plate scale is 5d = 1.5 mm per arcsec (R Friberg, 
private communication) so the illumination pattern on the primary 
moves many cm during routine operation, altering the fraction of 
the beam spilling off the edge of the mirror and terminating in the 
telescope enclosure. If dB/dR is the fraction of the total SCUBA 
response in an annulus of width dR at the edge of the primary of 
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diameter D = 15 m, the anticipated signal from variations in beam 
spill is 



dB ivDSd dB t , n -idB _ ± 
—— — — — — — — = 4 x 10 — — arcsec 

B tyD 2 /4 dR dR 



(6) 



This requires dB /dR ~ 1 per cent per mm, which is possible but 
larger than we expect. 

The secondary mirror support structure intercepts 8 per cent 
of the beam (P. Friberg, private communication). This coefficient 
could easily vary by one part in 10 4 per arcsec, which would pro- 
duce a signal at the observed level. We do not know whether this ef- 
fect produces a maximum or a minimum near the optic axis. In any 
case, this very large optical signal is not a malfunction of the tele- 
scope, but an unavoidable consequence of employing an observing 
strategy in which the se condary mirror is moved. 

iBorvs et"ai] l2004ll show that there is a spurious signal at a 
frequency near to 1 cycle per 16 samples present in SCUBA data 
collected during much of 2003, and other groups agree (A. Mortrier 
in preparation, T. Webb in preparation). Published power spectra of 
this spurious signal look strikingly similar to Fig.|5|given that the 
authors have not re-ordered their data to be in time-order before im- 
plementing a Fourier transformation. We presume that either hard- 
ware or software difficulties during that period lead to imperfect 
nod cancellation, leaving a substantial portion of the secondary 
mirror-induced pick up signal present in the nominally cleaned 
data. As a check, we have examined photometry data collected dur- 
ing the same period. In Photometry mode, SCUBA executes a small 
nine-point grid to account for small errors which might occur in the 
absolute pointing of the telescope. If the excess noise discussed in 
IBorvs etail l2004h is due to the secondary mirror and poor cancel- 
lation, there will be noise at one cycle per nine samples in photom- 
etry data. If the noise has a different origin it might well be at the 
same audio frequency in photometry mode as in jiggle mapping, 
i.e. one cycle per 16 samples. Fig.Qshows the power spectrum of 
a typical bolometer in the first set of photometry data from 2003 
we have chosen to examine, and the excess power is clearly near 
to the frequency of secondary mirror motion for photometry. This 
broad band of excess noise is absent in data from either 1998 or 
2000 which we examined. 

We conclude that the spurious signal from 2003 is in fact 
poor cancellation of secondary motion induced pick-up due to 
poor fine control of the secondary mirror's postion. This produces 
large residual signals in the double-differenced data streams. Be- 
cause the amplitude of the pick-up varies smoothly across the array 
(e.g. Fig.|4}, we have investigated the possibility of using the pat- 
tern as a template for removing the signal directly from the data 
streams. Unfortunately, no bolometer data seem to exhibit a pick- 
up signal completely described by a gradient term, and it appears 
that it is not possible to construct a time series that is purely due to 
the pick-up from this template. 

On a positive note, because this signal is bright and stable, 
it may be possible to use it to provide a short term gain moni- 
tor of very high signal to noise ratio in much the same way that 
Jarosik etal]<2003l) use the total power response of their radiome- 
ters to obtain a short term gain monitor for WMAP from the system 
temperature. However, a study of the stability of this signal and the 
reliability of its dependence on operating parameters is required to 
implement this strategy. 

We do not know how large this effect is 4 arcmin off the optic 
axis where SCUBA2 intends to have some sensitivity, but we pre- 
sume from Fig.[6jthat it is very large. Furthermore, it is unclear how 
well this signal can be disentangled from the atmospheric signal 
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Figure 7. The PSM of typical reduced photometry data taken between early 
to mid 2003 (bolometer 14 from run 57, 2003 February 9). The jiggle fre- 
quency spike is present (shown as a dashed line at / = 1/9 sample -1 ), as 
is a high level of broad-band noise between about 0.05 and 0.15 sample -1 . 
Both are caused by improper nod cancellation; the optical pick-ups are 
brighter in one nod position than the other, so the difference does not ef- 
fectively remove common mode contributions. While the photometry mode 
jiggle frequency is heavily contaminated, no spurious noise is present at 
the standard jiggle map frequency (the dotted line at / = 1/16 sample -1 ). 
This is further evidence that the pick-up is due to secondary mirror position. 



without a chopper present in the system. At this point, moving the 
secondary in any way seems to require nodding or some other sym- 
metric differencing technique to remove the pick-up signal. Unfor- 
tunately, the DREAM algorithm currently under consideration for 
use with SCUBA2 at the JCMT could well be affected by this sig- 
nal. Furthermore, instruments like THUM PER JWalker et alj2002l) 
and BOLOCAM iMauskop f et alj|200(j) . which may be deployed 
at the JCMT, will have to develop observation strategies which ac- 
count for this type of pick-up. These are challenging issues which 
must be addressed before the next generation of instruments is de- 
ployed at the JCMT. 
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